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Abstract 

We study the impact of the positivity of the measure of CP violation in neutrino 
oscillations on the values of \V e ^\ and |V^i|. We use the available information on 
the solar and atmospheric mixing parameters \V e 2\ and {Vu^l provided by neutrino 
oscillation data. We conclude that large values of \V e s\ are preferred at 95% confi- 
dence level. This, in turn, favours large values for Ji ep ton compared to CP violation 
in the quark sector. 



1 Introduction 



There are now compelling observations that neutrinos (leptons) do mix like quarks and 
have non zero masses 121 El S E] • Therefore we have a mixing matrix, analogous to that 
of the quark sector. In the absence of a fundamental theory of the origin of mass, this 
matrix is parametrised by several indenpendant parameters. Clearly, one of the important 
goals of particle physics in the next years will be to determine these parameters with the 
highest precision possible. One hopes that the knowledge of this matrix will point the 
way to a more fundamental theory of lepton mixings. Lepton mixing also has far reaching 
consequences in leptogenesis and in baryogenesis [7j. 

Over the last few years, several groups have strived to set bounds on the mixing 
parameters of the leptonic sector. Two preferred scenarios have emerged: the bi-maximal 
scenario [S] and the tri-bimaximal scenario 0. All these scenarios are summarized in 
reference jTU] . Recent data on neutrinos oscillations ^1] favour the tri-bimaximal scenario 
of lepton mixing. 

The purpose of this article is to study the role played by the measure of CP violation 
|12| . namely J, in the structure of the leptonic mixing matrix ^3]- A similar study has 
proven usefull before in the quark sector [Tlj. We will show that by using the positivity 
of J 2 , it is possible to restrict the range of | V^3 1 2 and |V^i| 2 in a usefull way. As we will 
see shortly, the maximum amount of CP violation in neutrino oscillations is gauged by 
the value of sin 6*13 [TH], which itself is closely related to |V e3 |. We will then concentrate on 
this last parameter and show that, at 95% confidence level (c.L), 0.0065 < | V^3 1 2 < 0.05. 

In the next section, we will describe the parametrization that we used for our analysis, 
which we will explain in the third section. In the fourth and fifth sections we will present 
our results on | V^=3 1 2 and |V^i| 2 , and on J 2 . We will conclude in the last section. 

2 Rephasing Invariant Parametrization of the Mix- 
ing Matrix 

We will assume that we have three neutrino flavors. In this case, the elements of the 
lepton mixing matrix are generally defined by 




(1) 



In terms of mixing angles, the mixing matrix becomes ^T] 
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where U is a Cabibbo-Kobayashi-Maskawa like unitary matrix and CP- violation a la Dirac 
is implemented through the phase 5. P is a matrix that contains Majorana CP violating 
phases. These phases are not relevant to our analysis as long as we consider only the data 
from oscillation experiments ^H]. Therefore, we can ignore them in what follows. 

In order to describe CP- violation in neutrino oscillations, the mixing matrix Vpmns 
can be parametrized in terms of four moduli ^7j; we choose the following four independant 
moduli: \V e2 \, \V e3 \, |V^i| and |V^ 3 |. 

As the Vij are complex quantities, setting V e3 = gives us two constraints: | V^3 1 = 
and 
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The whole mixing matrix is then described by \V e 2\ and 
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It is now possible to estimate the maximal value for \Vui\- In the bi-maximal scenario 
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In the tri-bimax scenario, one sets \ V f 
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We argue that this maximal value of |V^i| is quite reasonable as any non-zero value 
of K3 wm increase it or decrease it by an amount close to |K>3| or smaller IT7] . 
Nonetheless, we will study two cases later on: (V^il^, = 0.25 and 0.35. 
The exact expressions between the moduli of the elements of the mixing matrix (Vpmns) 



and the mixing angles experimentally accessible, namely 9 sun , 
in Xing [TO]; setting \ V e s\ = leads to IK2I = sin6 sun and [ 1 : 
that we need here are: 



6 at m and 6*13, can be found 
= sin9 atm . The expressions 
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sin 2 9i 3 
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Our current knowledge on the magnitude of the elements of the leptonic mixing matrix 
comes from experiments on neutrino oscillation and can be summarized as follows [THj: 



sin 2 0i 2 = 0.3 ±0.08 (9) 

sin 2 9 23 = 0.5 ±0.18 (10) 

sin 2 #i3 < 0.05 at 3a (11) 

Note that these values are close to those of the tri-bimaximal scenario. 
The measure of CP violation, J, is given in general as: 

j = %m(v e2 v, 3 v; 3 v; 2 ) (12) 



It has been shown before that if any element of the mixing matrix goes to zero, CP 
violation a la Dirac vanishes ^7j. However, the vanishing of CP violation does not imply 
that any element of the mixing matrix is zero. 
From the previous expression for J, we can infer an upper bound: 

\J\ < \V e2 \ \V^\ \V e3 \ [V^l (13) 



Using the values of the tri-bimaximal scenario (IK2I = l/v3 = (V^l, | 1 = 1/V2) we 
can get a rough estimate on the maximum value of |J|: 

\J\ < j$ (") 

Numerically, by taking the maximal value of jV^I ~ V0.05, we have | J\ ma x ~ 0.05 

As can be seen from eqs. 1 and 2, \V e3 \ has the simplest relation to Vpmns an d CP 
violation. It is also a measure of the maximum amount of CP violation possible in the 
leptonic sector. We will therefore pay special attention to this parameter. 
In terms of our 4 independant parameters, the measure of CP-violation is given by 



J 2 = (1 - |^ e3 | 2 - l^niv^l 2 !^! 2 ^! 2 
i (\v e2 \ 2 - + |v^| 2 |Ve3| 2 - iv^sl 2 !^! 2 - |v; 3 | 2 |v e2 | 2 ) 2 (15) 
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3 Analysis 



For our analysis, we use only the known mixing parameters: the solar mixing 612, the 
atmospheric mixing # 23 ■ We will use the positivity of J 2 to set an interval on | V^>3 1 2 and 
|V^i| 2 consistent with data. In the following, we will assume that cos 2 #i 3 ss 1 and later on 
we will estimate how a variation in this value could affect our results. This means that we 
use |V e2 | 2 ~ sin 2 9 12 = 0.3 ± 0.08 and |V^ 3 | 2 & sm 2 23 = 0.5 ± 0.18 and set the maximal 
value allowed for |V^i| at that obtained from the bi- maximal scenario: 

(16) 

We performed our numerical analysis as follows. Using the central values for |V e2 | 2 
and |V^ 3 | 2 we use eq.-15 and require that J 2 be positive. This gives us a range of values 
for \V e 3\ 2 and |V^i| 2 , as can be seen on fig. 1: any combination within the parabola leads 
to a positive value of J 2 . In order to take into account the uncertainties in \V e 2\ 2 and 
|V^ 3 | 2 , we used a Monte Carlo method. We generated values for \V e2 \ 2 and for |V^ 3 | 2 , 
independantly, according to the Gaussian distribution exp[— 0.5((x — x c )/<j x ) 2 ] centered 
at x c and with a standard deviation of a x . We covered the whole range from to 0.6 for 
1 1/^2 1 2 and to 1 for |V^ 3 | 2 and verified that our distributions were centered at 0.3 and 0.5 
with standard deviations of 0.08 and 0.18. For each combination of \V e 2\ 2 and |V^ 3 | 2 , we 
compute a curve similar to that of fig. 1. The numbers that we present here were obtained 
from 50,000 such combinations of (|V e2 | 2 , IV^] 2 ). Once we have these 50,000 curves, we 
scan the |V e3 | 2 — |V^i| 2 plane and count how many combinations of (|V e2 | 2 , (V^) 2 ) allow a 
given point on the grid; we used a grid of 40,000 points (200 on each axis). If a point is 
allowed by the lower branch of a parabola, it gains a weight of 1, similarly for the upper 
branch of a parabola. If a point of the grid is rejected by a lower branch of a parabola, 
it looses a weight of 1; similarly for the upper branch of a parabola. This results in a 
given point of the grid gaining a weight of 2 if it is allowed by a parabola and gaining a 
weight of if it is not allowed. This procedure leads to fig. 2-a where the probability of a 
given combination is given by the volume under the surface. The total volume under the 
surface is unity since we know that the final combination of parameters will be somewhere 
in our parameter space. The shape of this figure is easily understood: the base of the 
parabola will move a lot along the |V^i| 2 axis as values of \V e2 \ 2 and |V^ 3 | 2 are varied, 
thus excluding several points in the |V e 3| 2 _ l^i| 2 plane that are close to the (V^i) 2 axis, 
while points for larger values of |V e3 | 2 will be accepted by several parabolas. In order to 
put a confidence level on the combination (|V e3 | 2 , |V^i| 2 ), one simply slices the surface of 
fig. 2-a horizontally, starting from the top. Each slicing intersects the surface and the 
volume under the surface defined by this intersection represents the probability of finding 
the combinations (|V e3 | 2 , |V^i| 2 ). When this volume reaches 68% of the total volume, the 
curve that results of the intersection of the surface and the slicing is the 68% c.l. curve in 
the |V e3 | 2 and |V^i| 2 parameter space. The same applies for the 90% and 95% c.l. curves. 
This leads to figures 2-b. 
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In order to set c.l. limits on |V e3 | 2 alone, one has to integrate over |V^i| 2 . Starting 
with fig. 2-a, it is easy to integrate over |V^i| 2 : we simply add all the volume elements 
along an axis. This results in fig. 3, where the total surface under the curve is 1 since 
the physical, final parameter must be under this curve. In order to establish our c.l.s, we 
proceed as before and slice our figure horizontally, starting at the top. The intersection 
of our curve and the slicing, will define a certain range for |V e3 | 2 and when this range will 
cover a surface under the curve that represents 68% of the total surface, this range will 
be our 68% c.l. for |V e3 | 2 . The same procedure applies for the 90% and 95% confidence 
level ranges. 

4 Results 

According to this procedure, we can say that the values of \V e2 \ 2 = 0.3±0.08 and |V^ 3 | 2 = 
0.5 ±0.18 and the allowed ranges < \V e3 \ 2 < 0.05 and < \V^\ 2 < 0.25 lead to a 68%, a 
90%, and a 95% c.l. ranges in the | 1 2 — |V^i| 2 parameter space as shown on figure 2-b. 
Globally, |V^i| 2 has to be larger than 0.027 at 95% c.l. and \V e3 \ 2 has to be larger than 
0.0017, also at 95% c.l.. These bounds become 0.034 and 0.0033, respectively, at 90% c.l. 
If we intergrate over |V^i| 2 , we can say the 68% c.l. range for |V e3 | 2 is 0.0228 - 0.05; the 
90% c.l. range is 0.0103 - 0.05 and the 95% c.l. range is 0.0065 - 0.05. Therefore, we can 
say that values of | V^3 1 2 smaller than 0.0065 are excluded at the 95% c.l.. 

We can also project ourselves in the future and ask what these limits on | V^> 3 1 2 , after 
integration over |V^i| 2 , will become once the errors on \V e2 \ 2 and |V^ 3 | 2 will have been 
reduced. We performed the same analysis as above but we reduced the error on | V^2 1 2 
and I V^3 1 2 to 0.02 and 0.04, respectively. We studied 9 different combination of \V e 2\ 2 and 
|V m3 | 2 : (0.22, 0.30, 0.38) x (0.32, 0.50, 0.68), respectively. The bump that appears on fig. 
2-a moves substantially over the allowed ranges of |V e3 | 2 and |V^i| 2 but once we integrate 
over |V^i| 2 the allowed ranges of | V^> 3 1 2 change very little with the different combinations: 
0.023 - 0.05 at 68% c.l., 0.011 - 0.05 at 90% c.l. and 0.0068 - 0.05 at 95% c.l. The only 
exception being (0.38, 0.50) where they are 0.021 - 0.05, 0.0093 - 0.05, and 0.0058 - 0.05. 

One could ask how these bounds on \V e3 \ 2 will vary if we allow the ranges on \V e3 \ 2 and 
|V^i| 2 to vary. We explored this possibility and changed the ranges to < \V e3 \ 2 < 0.07 
and < |V^i| 2 < 0.35. From fig. 2-a, we can expect the general behaviour to continue 
for larger values of \V e3 \ 2 and |V^i| 2 . The different c.l. ranges obtained from the central 
experimental values of \V e 2\ 2 and |V^ 3 | 2 with their current errors, are shown on fig. 2-c. 
The different ranges are: 0.010 - 0.07 on \V e3 \ 2 and 0.051 - 0.285 on \V^i\ 2 at 68% c.l, 
0.0025 - 0.07 on |Ks| 2 and 0.021 - 0.333 on {V^ 2 at 90% c.l, 0.0011 - 0.07 on \V e3 \ 2 and 
0.014 - 0.35 on |V^i| 2 at 95%. Furthermore, when we integrate over (V^ij 2 , the c.l. ranges 
on |Ks| 2 are 0.0322 - 0.07, 0.0147 - 0.07, and 0.0091 - 0.07, at 68%, 90%, and 95%, 
respectively. This strenghtens our statement that values of \V e3 \ 2 smaller than 0.0065 are 
excluded at the 95% c.l.. 
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In this context of larger bounds on |V e 3| 2 and |V^i| 2 , we also projected ourselves in the 
future and redid the study with smaller errors on \V e2 \ 2 and |V^ 3 | 2 , namely 0.02 and 0.04. 
The results are that, again the different c.l. limits are rather insensitive to the different 
combinations of \V e 2\ 2 and |V^3| 2 used and read as follows: 0.033 - 0.07 at 68% c.l., 0.015 
- 0.07 at 90% c.l. and 0.0095 - 0.07 at 95% c.l. The only exception being (0.38, 0.32) 
where they are 0.030 - 0.07, 0.014 - 0.07, and 0.0088 - 0.07. 

5 What can we say about J 2 ? 

At this point, one can ask what the typical values of J 2 in the allowed ranges are. On figure 
4, we simply plot J 2 as a function of \V e ^\ 2 and |V^i| 2 for \V e2 \ 2 = 0.3 and |V^ 3 | 2 = 0.5; the 
intersection of the cone with the plane is our fig. 1. On this figure, the maximum value 
of J 2 is 2.42 x 10~ 3 . One also sees that J 2 decreases rapidly as we move away from the 
peak. Note however that for the allowed ranges discussed previously, the values taken by 
J 2 will be typically larger than 3 x 10 -4 . In order to give a better estimate, we generated 
combinations of our 4 parameters with the weights given by figure 2-a and calculated J 2 
for each combination. The results are that for all our c.l. ranges, J 2 takes values between 
~ 3 x 10~ 8 and 2.7 x 10~ 3 while the standard deviations of the distributions are constant 
at 0.6 x 10~ 3 . The average values of the distributions vary a little: 1.1 x 10 -3 over the 
68% confidence level range, 0.98 x 10~ 3 over the 90%, and 0.97 x 10~ 3 over the 95% 
cl range. Globally, these distributions do not vary much over our different c.l. ranges. 
Comparing figures 2-b-c with figure 1, it is clear that we cannot exclude J 2 = since the 
boundaries between J 2 positive and negative are well within our confidence level ranges 
for many combinations of \V e 2\ 2 and |V^3| 2 . Our real minimal value is in fact zero and 
only numerical accuracy made it non-zero. What is important here is not the minimal 
value but rather the average values over the different domains: these are clearly far from 
zero and we can say that smaller values of J 2 are disfavoured. 

As for J, the results are practically independant of the c.l. ranges and read as J m i n ~ 
2 x 10~ 4 , J max = 0.05, with an average value of 0.03 and a standard deviation of 0.01. 
The average values are the more meaningful results and small values are again strongly 
disfavoured. These are the strongest statements that we can make about J 2 and J. 

From the previous variations on our parameters, it is clear that our initial approxi- 
mation (ie \V e2 \ 2 ~ sin 2 6i2 = 0.3 ± 0.08 and |V^ 3 | 2 « sin 2 6 23 = 0.5 ± 0.18) has very little 
impact on our results. 

6 Conclusion 

The positivity of J 2 has proven a useful tool in setting constraints on different parameters 
of the leptonic mixing sector. Using only this variable, we can say that values of | V^3 1 2 
smaller than 0.0065 are excluded at 95%, stronly disfavoured, by current data on \V e2 \ 2 
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and | 1 2 . Such large values are in fact preferred by some authors as they can help in 
establishing the complementarity relations between quark and lepton mixing ^2]- From 
fig. 4, one sees that a large value of | V^3 1 2 is not a guarantee of CP violation but simply 
favours large CP violation. 

We can also say that J 2 must be smaller than 2.7 x 10~ 3 and that values much smaller 
than 1 x 10 -3 are strongly disfavoured. This translates into J ma x ~ 0.052 and values much 
smaller than 0.03 being also strongly disfavoured. Such values of Ji ep ton are typically 1000 
times larger than values of J quar k- 

When looking at figures 1 and 2, one can see that the positivity of J 2 will not be able 
to exclude | V^>3 1 2 = since the parabola will always extend all the way down to \V e z\ 2 = 0, 
no matter what the experimental precision on IK2I 2 and | 1 2 is. All we can do is set 
some confidence limits. 

Considering figures 1 and 4, it is clear that the tightest indirect contraints on | V^3 1 2 
and |V^i| 2 will come from an upper limit or a direct measurement of CP violation in 
neutrino oscillations; from a limit or a direct measurement of Jiepton- This would reduce 
the space to a narrow strip following the parabola of fig. 1. Short of such a measurement, 
one will have to rely on independant measurement of | V^3 1 2 and (V^il 2 . A measurement 
of |V^i| 2 far away from 0.15 (assuming that IK2I 2 and |V^3| 2 keep their current central 
values) would lead to the smallest possible range on | V^3 1 2 and push it far from 0, thereby 
increasing the probability of large CP violation. Similarly, a measurement of \V e z\ 2 with 
a very small value would put the tightest limits on |V^i| 2 and at the same time would 
reduce greatly the possibility of large value of \J\ 2 . Assuming that a small value of | V^ 3 1 2 
will be difficult to measure, one concludes that experimental efforts should be devoted to 
measuring |V^i| 2 since this could potentially strongly restrict the range of | V^>3 1 2 . Let's not 
forget that complete determination of the neutrino mixing matrix will shed some light on 
the deep structure of the leptonic spectrum and help us understand better the properties 
of neutrinos and mass generation in general. 
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8 Figures Captions 

Figure 1: Limits in the \V e3 \ 2 — |V^i| 2 plane obtained from the positivity of J 2 . The 
parabola represents J 2 = 0; the space inside represents positive values of J 2 . 
Figure 2-a: Probability distribution in the | V e3 \ 2 — \ V^i \ 2 plane obtained from 50,000 curves 
similar to those of figure 1. The vertical axis is in arbitrary units as the volume under 
the curve, representing the total probability, is unity. 

Figure 2-b: Confidence levels obtained from figure 2-a. The 68% c.l. has the smallest 
surface while the 95% c.l. has the largest. 

Figure 2-c: Similar to figure 2-b but with < \V e3 \ 2 < 0.07 and < |V^i| 2 < 0.35. 
Figure 3: Probability distribution of | V^ 3 1 2 after integrating the [V^iJ 2 axis of figure 2-a. 
Figure 4: J 2 for different values of | V^ 3 1 2 and [V^i) 2 using [V^l 2 = 0.3 and [V^) 2 = 0.5 
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